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Abstract. We report on work done by the Virgo consortium, an interna- 
tional collaboration set up in order to study the formation and evolution 
of Large Scale Structure using N-body simulations on the latest generation 
of parallel supercomputers. We show results of 256 3 particle simulations of 
the formation of clusters in four Dark Matter models with different cos- 
mological parameters. Normalizing the models such that one obtains the 
correct abundance of rich clusters yields an interesting result: The peculiar 
velocities of the clusters are almost independent of O, and depend only 
weakly on T, the shape parameter of the power spectrum. Thus, it is nearly 
impossible to distinguish between high and low Q models on the basis of 
the peculiar velocities. 
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1. Introduction 

Clusters of galaxies are the largest objects one finds today in the Universe. 
They are well studied from both observational and theoretical viewpoints 
(for a recent review c.f. Bahcall (1996a)). 

It may be possibile to distinguish between high and low 9 Universes 
by using the peculiar velocities of the clusters. On this basis Bahcall et 
al. (1996b) concluded that the observed clusters velocity function "is most 
consistent with a low mass-density ($1 ~ 0.3) CDM model". However, we 
show in section 4 that this is dubious when models are normalized such 
that they all give the correct abundance of rich clusters. We present results 
from the simulations that support our theoretical considerations in section 
5. 

Before we talk about the details of the analysis of clusters from the 
simulations we give in section 2 a brief description of the numerical code. 
In section 3 we discuss the models used. 

2. The Code 

The Virgo consortium was formed in order to study the evolution of both 
Dark Matter and gas in the expanding Universe using the latest generation 
of parallel supercomputers. 

The code used is called "Hydra" which was developed by Couchman et 
al. (1995) and parallelized by Pearce et al. (1995). Gravity is treated using 
an adaptive P 3 M technique. This places mesh refinements on the regions 
of strongest clustering. Large refinements are done in parallel across the 
processors, smaller ones are farmed out to single processors. 

The simulations were run on the Cray T3D supercomputers at the com- 
puter center of the Max Planck Society in Garching and at the EPCC in 
Edinburgh. 

3. The Dark Matter Simulations 

We have carried out a set of very large N-body simulations of CDM uni- 
verses with four different choices of parameters. The models chosen were 
Standard CDM (SCDM), a high Q model with an additional radiative com- 
ponent (rCDM), an Open CDM model (OCDM) and a flat low 9. model 
with a Cosmological Constant (ACDM). With the exception of the Open 
Model, which only had 200 3 particles, each simulation followed the evolu- 
tion of 256 3 particles in a box of 239.5 h~ 1 Mpc 1 on a side. 

As usual we express the Hubble constant as Hq = 100 h km/ (Mpc sec). 



The Virgo Consortium: The evolution and formation of galaxy clusters 
TABLE 1. The Virgo models 



Model fl A h a s 



SCDM 


1.0 


0.0 


0.5 


0.51 


0.50 


rCDM 


1.0 


0.0 


0.5 


0.51 


0.21 


ACDM 


0.3 


0.7 


0.7 


0.90 


0.21 


OCDM 


0.3 


0.0 


0.7 


0.85 


0.21 



In all models, the initial fluctuation amplitude was set by requiring that 
the models should reproduce the observed abundance of rich clusters (for 
details see section 4). Table 1 gives the parameters of the models. 

Figure 3 shows the evolution of the same cluster 2 in the rCDM and 
the ACDM model for the redshifts 2, 1, and 0. Structure forms earlier in 
the low Q universe. But already at a redshift of 2 a very large filamentary 
structure can be seen in the high $7 universe, too 3 . 

4. Peculiar Velocities of Galaxy Clusters from Linear Theory 

According to linear theory, the mean-square peculiar velocity of galaxy 
clusters is 

1 f 00 

{v * D) = 2^ H ° nl 'Jo P ^ T ) w2 ^ R ) dk > C 1 ) 

where W 2 (k, R) = exp[— k 2 R 2 ] is a window function of radius R (see, 
e.g., Croft & Efstathiou (1994)). The power spectrum P(k) is taken in 
its parametric form introduced by Bond & Efstathiou (1984), 

P ^ = {1 + [ak + {bkfl 2 + (cAi) 2 ]"} 2 /" ' ^ 

where a = (6.4/T) h' 1 Mpc, b = (3.0/T) h' 1 Mpc, c = (1.7/T) h' 1 Mpc and 
v = 1.13. The quantity V is given by 

_ J S\)/i models with Oo + Ao = 1 ,„, 

" | f)oV[0-861 + S^mioTd) 2 / 3 ] 1 / 2 models with decaying v '' 

where mio is the neutrino mass in units of lOkeV and Td is its lifetime in 
years (Bond & Efstathiou (1991)). Equation (3) is valid for any model in a 
spatially flat universe. 

The simulations were run with the same phases. 
3 The clusters were taken from an additional set of runs with the same parameters 
as above but box sizes of 85 h^ 1 Mpc and 141 h Mpc for the high and low £1 models, 
respectively. 
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Figure 1. The evolution of the same cluster in the rCDM (leftmost pictures) and in the 
ACDM model (rightmost pictures) , for redshifts of 2 (top) , 1 (middle) , and (bottom) . 
The sizes of the regions shown are 21 x 21 x 8h _1 Mpc 3 and 35 x 35 x 13 ft -1 Mpc 3 for 
rCDM and ACDM, respectively. 
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The normalisation B of the power spectrum can be obtained in two 
different ways: Either by using a relationship between B and the COBE 
measurements (as shown in Efstathiou et al. (1992)) or by relating B to 
the rms linear fluctuation in the mass distribution on scales of 8/i _1 Mpc, 
erg, which is defined by 



where R% = 8 /i _1 Mpc, and ji is a spherical Bessel function. 

Values for og were obtained by either using the mass or the X-ray 
temperature functions of rich clusters (White et al. (1993) (WEF), Viana 
& Liddle (1996), Eke et al. (1996) (ECF)). WEF obtain a 8 « 0.57ft/- 56 , 
ECF get 

_ / (0.52 ± 0.04) Q ^e+o.wn for Aq = Q 
a8 ~\ (0.52 ±0.04) ft-°- 52+a13 a ° for ft + A = 1 ( j 

with the quoted statistical uncertainties obtained using bootstrap methods. 
Inserting eq. (4) into eq. (1) yields 

{vl D ) = W 4 (n°- 6 a 8 ) 2 f(T,R) [km/sec], (6) 

with the abbreviation 

HT.R )S ^ SnW'HKRyik (7) 

/f(*)(jfeji (*«»)) <Pk 

Comparing eqs. (5) and (6) one sees that (ft a6 <78) 2 only very weakly de- 
pends on ft. Thus, the peculiar velocities of galaxy clusters are nearly a 
function of T alone. 

Taking R = 1.5 hr 1 Mpc one obtains 

/(r, R = 1.5 /i" 1 Mpc) = 12.5 r~ L08 + 49.4 (8) 



5. Peculiar Velocities of Galaxy Clusters from the Simulations 

When comparing N-body simulations with real data one has to find a 
way to select the clusters in the simulation. WEF find clusters by locat- 
ing high-density regions with a friends-of-friends group finder, and then 
getting masses from the particle count within spheres of comoving radius 
r = 1.5 h^ 1 Mpc. We use the same algorithm and treat objects with a mass 
larger than 5.5 • 10 14 /i _1 Mq as galaxy clusters. 
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TABLE 2. Peculiar velocities of galaxy clusters 
Model {vt D )u [km/sec] {vl D )sim [km/sec] JV 



SCDM 


441 


461 


20 


rCDM 


540 


492 


21 


ACDM 


462 


452 


21 


OCDM 


437 


476 


18 



In table 2 we give the peculiar velocities for the four models from linear 
theory (first column, using eq.s (6) and (8)), from the simulations (second 
column), and the number of clusters found in the simulations (third col- 
umn). Linear theory clearly predicts the simulation results to within their 
uncertainties. 

From the numbers it is obvious that the peculiar velocities are not a 
good way to discriminate between the models if these are normalized such 
that they reproduce the correct abundance of rich clusters. In particular, 
there is no difference between high and low f2 models. 
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